Bound spinons in an antiferromagnetic S=| chain with a staggered field 



M. Kenzelmann 1 ' 2 , Y. Chen 1 ' 3 , C. Broholm 1 * 2 , D. H. Reich, 1 and Y. Qiu 2 < 4 
(1) Department of Physics and Astronomy, Johns Hopkins University, Baltimore, Maryland 21218 
(2) NIST Center for Neutron Research, National Institute of Standards and Technology, Gaithersburg, Maryland 20899 

(3)Los Alamos National Laboratory, Los Alamos, NM 87545 
(4)Department of Materials Science and Engineering, 
University of Maryland, College Park, Maryland, 20742. 
(Dated: February 2, 2008) 

Inelastic neutron scattering was used to measure the magnetic field dependence of spin excitations 
in the antiferromagnetic S=| chain CuCl2 • 2(dimethylsulfoxide) (CDC) in the presence of uniform 
and staggered fields. Dispersive bound states emerge from a zero-field two-spinon continuum with 
different finite energy minima at wave numbers q—n and q t ~ n(l — 2(S Z )). The ratios of the field 
dependent excitation energies are in excellent agreement with predictions for breather and soliton 
solutions to the quantum sine-Gordon model, the proposed low-energy theory for S=^ chains in 
a staggered field. The data are also consistent with the predicted soliton and n=l,2 breather 
polarizations and scattering cross sections. 
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Shortly after the advent of quantum mechanics, Hans 
Bcthc introduced a model antifcrromagnet that continues 
to play a central role in quantum many body physics . 
The isotropic antiferromagnetic (AF) spin-1/2 chain, has 
a simple spin Hamiltonian: 7i = J J2 n • S n +i, is inte- 
grable through Bethe's Ansatz, and is realized with high 
fidelity in a number of magnetically anisotropic Cu 2+ 
based materials. Because it sits at the boundary between 
quantum order and spin order at T = 0, Bethe's model 
is ideally suited for exploring quantum critical phenom- 
ena and the qualitatively different phases that border the 
critical point 0, Q • This letter presents an experimen- 
tal study of the profound effects of a symmetry breaking 
staggered field on excitations in the spin-1/2 chain. 

In zero field, the fundamental excitations of the spin- 
1/2 chain are not spin waves but domain wall like quasi- 
particles called spinons that separate reversed AF do- 
mains 0, HI El • The ground state is a Luttinger liquid 
and the spinons are non-local spin-1/2 objects with short 
range interactions. Thus, spinons can only be excited in 
pairs and produce a gapless continuum. Such a spectrum 
has been observed in several quasi-one-dimensional spin- 
1/2 chain systems @,@>I3 and is now understood to be 
a distinguishing attribute of quantum-critical systems. 
The dramatic effect of a staggered field was discovered 
through a high field neutron scattering experiment on 
the quasi-one-dimensional spin-1/2 antiferromagnet cop- 
per benzoate [l(i| . Designed to verify theoretical predic- 
tions of a field driven gapless incommensurate mode 
this experiment instead revealed a field induced gap in 
the excitation spectrum. The critical exponent of ~ 2/3 
describing the field dependence of the gap in copper ben- 
zoate, [PMCu(N0 3 ) 2 • (H 2 0) 2 ] n HU, and Yb 4 As 3 0, 
identified the source of this gap as the staggered field 
that accompanies a uniform field in materials with alter- 
nating Cu-coordination. The staggered field yields an 
energetic distinction between reversed domains, which 



confines spinons in multi-particle bound states [13j. 

A quantitative theory for this effect was developed by 
Affleck and Oshikawa starting from the following exten- 
sion of Bethe's model Hamiltonian jl^L fT3| . 

n = Si s i+1 + ^(-i) j 'd • (Sj_i x s 3 ) 

i 3 

-E H a [9^ + (-iyg s aff ]SP. (1) 

The alternating spin environment is represented by the 
staggered Dzyaloshinskii-Moriya (DM) interaction and 
Zeeman terms. Through an alternating coordinate trans- 
formation, the model can be mapped to a spin-1/2 chain 
in a transverse staggered field that is proportional to the 
uniform field H . While the zero field properties of Eq.(l) 
are indistinguishable from Bethe's model, an applied field 
induces transverse AF Ising spin order and a gap. Using 
bosonization techniques to represent the low energy spin 
degrees of freedom, Affleck and Oshikawa showed that 
their dynamics is governed by the quantum sine-Gordon 
model (QSG) with Lagrangian density 

£ = \[{dtW - (dxW] + hCct*(fl$). (2) 

Here h oc H is the effective staggered field. C(H) and 
(3(H) (which goes to \/2w for H — > 0) vary smoothly 
with the applied field and can be determined numerically 
through the Bethe Ansatz for h « H 

With applications from classical to particle physics, 
the SG model plays an important role in the theory of 
non-linear dynamic systems [lfij . Excitations are com- 
posed of topological objects called solitons that encom- 
pass a localized ±2ir/0 shift in <p for a soliton and anti- 
soliton respectively [17| . In addition there are soliton- 
antisoliton bound states called breathers, which drop be- 
low the soliton-anti-soliton continuum as the non-linear 
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FIG. 1: Normalized scattering intensity I(q,oj) in zero field 
and H—ll T in the S = ± chain CDC as a function of the chain 
wave- vector transfer q and energy transfer Tiw, measured at 
T=40 mK. The data were binned and smoothed, leading 
to an effective wave-vector resolution 5g=0.0577r and 0.061-7T 
and an energy resolution 0.1 and 0.18 meV for configurations 
A and B, respectively. The colorbar indicates the scattering 
strength. 

term is increased. The excited state wave functions are 
known for both solitons and breathers, which enables ex- 
act calculation of the inelastic scattering cross sections. 
In this Letter we use neutron scattering from a magne- 
tized spin-1/2 chain with two spins per unit cell to test 
these results and more generally to explore the dynamics 
of spinons with long range interactions. 

Based on the temperature dependence of the suscepti- 
bility and specific heat, CuCl 2 • 2((CD 3 ) 2 SO) (CDC) was 
identified as an AF S=\ chain system with J=1.5 meV, a 
staggered (/-tensor and/or DM interactions [lillllSj. The 
spin chains run along the a-axis of the orthorhombic crys- 
tal structure (Prima) [2(J, with the Cu 2+ ions separated 
by 0.5a ± 0.22c. Wave vector transfer is indexed in the 
corresponding reciprocal lattice Q(hkl) = ha*+kb*+lc* , 
and we define q = Q • a. Due to weak inter-chain inter- 
actions, CDC has long-range AF order in zero field be- 
low T/v=0.93 K with an AF wave- vector Q m = a*. An 
applied field strongly suppresses the ordered phase [l9| . 
indicating that inter-chain interactions favor correlations 
that are incompatible with the field-induced staggered 
magnetization [2 1| . Above the H c = 3.9 T critical field 
for Neel order, we find that CDC is an excellent model 
system for our purpose. 



Deuterated single crystals were grown through slow 
cooling of saturated methanol solutions of anhydrous 
copper chloride and deuterated dimethyl sulfoxide 
((CD 3 ) 2 SO) in a 1:2 molar ratio 0. The sample stud- 
ied consisted of four crystals with a total mass 7.76 g. 
The experiments were performed using the disk chopper 
time-of-flight spectrometer (DCS) at the NIST Center 
for Neutron Research with the c-axis and the magnetic 
field vertical. In configuration A the incident energy was 
i?i=3.03 meV and the a-axis was parallel to the inci- 
dent neutron beam direction k^. Configuration B had 
i?i=4.64meV and Z(ki,a)=60°. The counting time was 
18 hrs at 11 T and an average of 5 hrs for each measure- 
ment between and 8T. The raw scattering data were 
corrected for a time-independent background measured 
at negative energy transfer, for monitor efficiency, and 
for the Cu 2+ magnetic form factor, folded into the first 
Brillouin zone, and put onto an absolute scale using the 
elastic incoherent scattering from CDC. For the normal- 
ization, the H/D ratio (=0.02) was measured indepen- 
dently through prompt-7 neutron activation analysis. 

Figures ^a) and QJb) show that for T <C J/k Bl the 
zero-field excitation spectrum of CDC consists of contin- 
uum scattering above a low-energy threshold that varies 
as Ticu—^Jl sin(g)| through the zone 0. An exact analyti- 
cal expression for the two spinon contribution to the scat- 
tering cross section which accounts for 72.89% of the total 
spectral weight was recently obtained by Bougourzi et al. 
[ol I22I Ei^ . Figures El and 01 show a quantitative compar- 
ison of this result (blue line), duly convoluted with the 
experimental resolution, to the experimental data. The 
excellent quantitative agreement between model and data 
provides compelling evidence for spinons in the zero field 
state of CDC. Note that the Goldstone modes that are 
expected due to Neel order for hio < k^T^ ~ 0.1 meV, 
are not resolved in this experiment. 

Figures njc) and^d) show that the magnetic excita- 
tions in CDC change dramatically with field and are dom- 
inated by resolution-limited modes for H = 11 T. Figures 
EJandEJb) show spectra at the wave vectors correspond- 
ing to the minima in the dispersion relations, which oc- 
cur at q = 7r, and qi = 0.777T, as determined from the 
constant-?iw cut in Fig. Eta) . These data graphically il- 
lustrate the field-induced transfer of spectral weight from 
the two-spinon continuum into single-particle excitations. 
A phenomenological cross-section of long-lived dispersive 
excitations was fit to the data near q — tt and qi to take 
into account the experimental resolution and thereby ac- 
curately locate the excitation energies. These fits are 
shown as red lines in Figs. El and El and the inferred pa- 
rameters characterizing the dispersion relations are dis- 
played for a series of fields in Fig. El 

We now examine whether our high field observations 
are consistent with the QSG model for spin-1/2 chains 
in a staggered field 14]. First, the model predicts single 
soliton excitations at <7i=7r(l — 2(S Z )). This is qualita- 
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FIG. 2: Spectra at q—TT observed with configurations A (a) 
and B (b) for zero field and 11 T obtained from data in the 
range 0.957T < q < 1.057T. The zero-field intensity for 0.37T < 
q < 0.4-7T was subtracted as a background. The top axes 
indicates wave-vector transfer, k, perpendicular to the chain. 
The blue line is a fit of the two-spinon cross-section to the 
zero field scattering [(J- The red lines are fits to the H — 11 
T data as explained in the text, including the theoretically 
calculated breather continua polarized perpendicular to the 
field |15|. The inset in (b) highlights the high energy range. 

tively consistent with the raw data in Figs. 1(c) and 1(d), 
and with Fig. 4(d), which shows how qi moves across the 
zone with H. Quantitative agreement is also apparent 
from the solid line in Figure 4(d), which is the predicted 
field dependence as calculated from the magnetization 
curve for a spin- 1/2 chain @. 

The soliton and antisoliton mass is related to the ex- 
change interaction of the original spin chain and the uni- 
form and staggered fields, H and h, as follows 0, 0] 



M » J(2HB^) (1+0/2 x 



(3) 



Here £ = /3 2 /(8tt - /3 2 ) 1/3 for H -> and B = 
0.422169. Assuming h oc H, the soliton energy versus 
field is shown as a solid red line in Fig. EI a ) ■ While Eq. 
© is at the limit of validity for CDC at H = 11 T, we 
attribute the discrepancy with the mode energy at qi (red 
triangles) to inter-chain interactions that suppress the 
effective staggered field close to H c . The more general 
expression, h oc (H — H c ) a , yields a good fit for a — 
0.68(5) (dashed red line). 

The QSG model predicts that breather bound states 
of 2n-solitons should be accessible at q=n with masses 

M„ = 2Msin(n7r£/2). (4) 

Sharp modes are indeed observed in CDC at q = tt. Their 
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FIG. 3: Constant-oj (a) and constant-g scans (b) through 
the gapped soliton at H—ll T at the incommensurate wave- 
vector qi. The zero- field scattering is shown for comparison 
together with a fit of the two-spinon cross-section to the data. 
The unfitted peak corresponds to a moving breather. These 
data are from configuration B. The ranges of integration were 
(a) 0.675 meV < Twj < 0.75 meV, and (b) 0.72tt < q < 0.82tt. 



energies are compared to the breather masses predicted 
for h oc H (solid lines) and h oc (H - H c ) am (dashed 
lines) in Fig.^Ja). The ratios of the commensurate and 
putative breather mode energies to the lowest energy in- 
commensurate and putative soliton mode energy, shown 
in Fig.^Jb), are in excellent agreement with the normal- 
ized breather energies M n /M for n—1 and n=2. This 
comparison, which is insensitive to the origin of the stag- 
gered field, suggests that breathers indeed exist in CDC. 

The evidence for breathers is strengthened as we exam- 
ine the polarization of the scattering at q = ir. According 
to the QSG model, n=odd (even) breathers are polarized 
in the plane normal to H and perpendicular (parallel) to 
h|T^. Neutron scattering probes the projection of spin 
fluctuations on the plane normal to the scattering vec- 
tor Q. Figure [21 shows that the Tilo—1 meV peak seen 
for Qa ~ (1,1.64,0) in configuration A is absent for 
huj=l meV and Qb ~ (1,0,0) in configuration B. In 
a quasi-one-dimensional system the only explanation for 
this is that the excitation is polarized along Qb | |a| |h as 
expected for an even numbered breather, and hence is ex- 
tinguished by the polarization factor in the neutron scat- 
tering cross section for configuration B. The predicted b 
and c axis polarizations respectively of the n—l breather 
and the soliton are confirmed by the consistent polar- 
ization factor corrected intensities from configurations A 
and B for H = 11 T in Fig. He). 

One of the unusual aspects of the QSG model is that 
complex features such as the breather and soliton struc- 
ture factors can be calculated exactly [l5|. The solid 
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FIG. 4: (a) Field dependence of the soliton and breather 
modes. The lines are the QSG predictions for h oc H (solid) 
and h=c(H - H c ) (>8 (dashed), (b) Ratio of breather energies 
to soliton energies versus H, compared to the QSG model, (c) 
Energy integrated intensities divided by the QSG predicted 
polarization factors and compared to structure factor predic- 
tions using particle energies given by the dashed line in (a) 
and a common overall scale factor 15]. (d) Incommensura- 
tion qi obtained from the minimum in the soliton dispersion 
compared to the predicted qi — ir(l — 2{S z )) using the magne- 
tization curve (Sz) = l/w arcsin((l — n/2 + it J / (g z H ^b))" 1 ) 
for a Heisenberg spin-1/2 chain|5|. The field dependence was 
measured with configuration A, and an additional measure- 
ment at H = 11 T was performed using configuration B. 

lines in Fig.^c) show that these exact results are consis- 
tent with the field dependent intensities of the commen- 
surate and incommensurate low energy modes in CDC. 
For H = 11 T the third breather is expected at about 
1.4 meV, close to the energy of the soliton mode at h=l. 
The peak close to 1.4 meV has intensity / = 0.14(3) 
in configuration A and / = 0.26(2) in configuration B 
and this is consistent with a third breather contribution 
polarized along b. The inferred b-polarizcd intensity of 
/° xp — 0.23(7) is however much greater than the inten- 
sity predicted for the n = 3 breather (i^ SG = 0.026(7)), 
which indicates additional anisotropic contributions to 
the inelastic scattering there. 

In addition to the field-induced resonant modes, the 
experiment shows that a high energy continuum persists 
for ff=ll T. Fig.^d) for example clearly shows a broad 
maximum in the q— dependence of neutron scattering for 
energies hu> > 1.6 meV and q ~ n. Firm evidence for 
continuum scattering comes from the field dependence of 
the first moment (huj)Q=h 2 J 5*(Q, u))uidu), which is pro- 
portional to the ground state energy (Ti.) with a negative 
(/-dependent prefactor [24| . At zero field the experimental 
value of {Tilo)q corresponds to (H) = —0.4(1) J, in agree- 
ment with Bethe's result of (H) = (±-ln2)J rs -0.44J. 
At 11 T, however, (Ti) derived solely from the resonant 
modes is —0.25(6) J when (H) is expected to be — 0.34J 
|lf)j . The discrepancy is an independent indication of 
spectral weight beyond the resonant modes. For q = n, 
the transverse contribution to the continuum scattering 



predicted by the QSG model [l5j and shown as a solid 
line in the inset of Fig. [3 has a maximum close to a weak 
peak in the measured scattering intensity. The shortfall 
of the theoretical result suggests that there are additional 
longitudinal contributions to the continuum scattering. 

In summary, staggered field induced spinon binding in 
spin-1/2 chains provides an experimental window on the 
unique non-linear dynamics of the quantum sine-Gordon 
model. Our neutron scattering experiment on quasi-one- 
dimensional CDC in a high magnetic field yields clear 
evidence for soliton/ antisoliton creation at wave vector 
transfer q = ir(l — 2(S Z )), as well as n = 1 and n = 2 
breather bound states at q — tt. Interpretation of the 
data throughout the Brillouin zone will require exact di- 
agonalization studies and a better understanding of lat- 
tice effects than provided by the continuum field theory 
reviewed in this paper. Other results that call for further 
experimental and theoretical work are the observation of 
high energy continuum scattering in the gapped phase 
and indications that inter-chain interactions can renor- 
malize the soliton mass. 
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